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Abstract 

We investigate the implications of the recent measurement of muon anomalous magnetic 
moment for the direct detection of neutralino dark matter in the three different SUSY 
models: mSUGRA, a model with non-universal Higgs mass, and an 50(10) GUT model. 
We consider two extreme scenario for Aa M bound, i.e. 27 x 1(T 10 < Aa M < 59 x 1(T 10 
(la bound) and < Aa^ < 11 x 10~ 10 (2a below). In mSUGRA model, the counting 
ratio may be above the sensitivity of the future experiments when parameters are within 
1(7 bound of Aa^. However, the H x tends to be high compared to the currently favored 
value f2 = 0.3. For models with the non-universal scalar masses, the possibility to have 
the consistent £l x and the high counting ratio is open up in the region of parameter 
space where Higgsino mass \x is smaller than mSUGRA prediction. In particular, in 
the 50(10) model, the LSP dark matter detection rate may be enhanced by almost 
one order of magnitude compared to mSUGRA and the model with non-universal Higgs 
mass, for cosmologically acceptable VL x h 2 . The highest detection rate of LSP dark matter 
occurs in the region where the LSP constitutes a subdominant part of local halo DM. 
Implication of SUSY mass parameter measurement under the cosmological constraint is 
also discussed. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) [I]] is one of the best motivated exten- 
sions of the Standard Model. It offers a natural solution of the hierarchy problem M as well as 
amazing gauge coupling unification || . Since naturalness requires that at least some superpar- 
ticles have masses at or below the TeV scale, supersymmetric theories generally predict a rich 
phenomenology at future colliders such as Tevatron, LHC, or proposed linear colliders (LC) 
[f|]. As an extra "bonus", the simplest version of the MSSM, where R— parity is conserved, 
also contains a new stable particle (the lightest supersymmetric particle, LSP); in most cases 
this is the lightest neutralino, which often makes a good Dark Matter(DM) candidate ||. 

Though we might have to wait for LHC experiment to start for the direct discovery of 
supersymmetry, we may be able to probe SUSY models earlier by using precision measurements 
in low energy experiments. Sparticles contribute to low energy measurements through loop 
effects, and these effects may become significant if the their masses are not too large. Another 
possibility is direct detection of dark matter by the detector placed deep underground || [7|, 
||, [T(J . Note that the existence of dark matter already indicates the needs of new physics. 

Recently the Brookhaven E821 experiment has released a measurement of the muon anoma- 



lous magnetic moment, reporting a 2.6a deviation from the standard model value |IT] 

Aa M = a e * p - af* = 43(16) x 1(T 10 , (1) 

which has generated considerable interest. This deviation may indeed be a sign of new physics 
beyond the standard model (SM) and could be accommodated by supersymmetric contribu- 
tions |I~2]l- |f26[1 . There are essentially two types of diagrams involving superparticles which 
contribute to a p , i.e., neutralino-smuon and chargino-sneutrino loop diagrams |27| . Since a p 
requires chirality violation, for tan/5 ^> 1 the dominant contributions are proportional to the 
product of an electroweak gauge coupling and the Yukawa coupling of the muon, where the 
latter factor either comes directly from the higgsino component of the chargino or neutralino 
in the loop, or from JIl—P-r mixing. For tan/3 3> 1 one has thus Aa M oc tan/3. If mass splittings 
between different sparticles are not too large, so that their mass scale can be described by the 
single parameter m SU sY, the overall result can be estimated as [^| 

where g\ and g 2 are the U{l)y and SU{2) gauge couplings, respectively. In our convention, the 
sign of Aa^ is equal to the sign of Higgsino mass parameters fi, when the gaugino masses and 
tan/3 are positive. Therefore Aa^ imposes non-trivial constraints on the MSSM parameters 
such as the sign of fi, superparticle masses and the ratio tan/3 of the vacuum expectation 
values (vevs) of the two neutral Higgs fields of the MSSM. The new data prefers a positive /i 
and one has an upper or lower limit on the chargino and sneutrino masses for a given tan/3. 
Through these non-trivial constraints, Aa^ will affect the predictions of other observables. 

There exist relations between the SUSY contribution to a p and the predicted LSP-nucleon 
scattering cross section |2£J. A direct connection comes from the common dependence on 
tan/3, since both a M and o XP increase with increasing tan/3. Similar to the SUSY contribution 
to a M , the leading contributions to the spin-independent (coherent) contribution to a xp also 
involve violation of chirality [B0|. Main contributions come from the exchange of CP-even 



Higgs bosons. It couples to the strangeness component of proton or neutron, whose coupling 
has the same tan/3 dependence as the Yukawa coupling of the muon. Note also that the sign 
of /i, favored by measurement, suggests large value of a XP for smaller LSP mass due to the 
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lack of accidental cancellation of couplings [ 29| . Importance of DM search experiments have 
now increased substantially. 

In general, larger o XP implies smaller relic DM density Q x h 2 . Since the Higgs-LSP-LSP 
couplings require higgsino-gaugino mixing, they scale like for /i 2 3> M|. The Higgs 
exchange contribution to a xp then scales like l//i 2 . On the other hand, the change in the 
size of \x also affects Vt x h 2 through the LSP dark matter annihilation into gauge bosons or 
s-channel Higgs boson. The annihilation cross section for XX ~* WW and ZZ are highly 
suppressed in the bino-like LSP case which is favored in mSUGRA model. However, a large 
cross section is possible when the LSP has a significant Higgsino component. This is because 
WW cross section depends on Wxx ± interactions and only couplings WH°H ± and WW°W ± 
are allowed by gauge invariance. For ZZ cross section, it depends on Z\X°i interactions which 
is possible only through ZH°H° couplings. Hence a reduction of fi (an increase in higgsino 
component of the LSP) implies an increase of LSP annihilation cross section, which in turn 
gives a reduction of Q x h 2 . The pair annihilation into fermions through s-channel pseudoscalar 
Higgs boson is also enhanced if LSP is the mixture of gaugino and higgsino. 

The relations between and a xp , and a xp and Q, x h 2 we discussed so far are model inde- 
pendent. One may further investigates more model dependent aspect of the relation. This has 
been intensively studied in a class of model called mSUGRA with an universal scalar mass 
m, an universal gaugino mass M, and an universal trilinear coupling A |29[ 20, 21, |26[. The 



model relates a lower bound of higgs mass and Aa^ to the constraints to right handed slepton 
and LSP masses, and a simple flavor structure of soft masses leads the prediction to b — > 57. 

The parameter space of mSUGRA is then quite limited. This is because the model predicts 
\i 3> Mi unless tan j3 ^> 1 or m ^> M . This tends to lead to too large a fl x h 2 , because the LSP 
pair annihilation cross section, a xp (xx ~* X) is proportional to m x /m~ while Higgs mass and 
a M constraints give a lower bound of those particles masses in general. The b — > S7 constraint 
also pushes up the overall SUSY scale. Especially for large tan/3 (~ 30), one might need 
stau-neutralino coannihilation or a s-channel annihilation of the neutralino to make relic DM 
density reasonable ||31|| . The DM signal also turns out to be rather low; the maximal value of 
the cross section is only slightly above the proposed sensitivities. 

Note however that the constraint strongly relies on the assumption that all scalar masses 
are universal at GUT scale. The rather parameter independent prediction /i > M in mSUGRA 
is actually relying on the strict universality. Smaller value of fi can be easily achieved in a 
model with non- universal Higgs masses | ]29[| . This leads to consistent DM density and larger 
counting ratios. The constraint from b — > 57 is sensitive to on the structure of scalar masses. 

A study of non-mSUGRA model may be useful in this stage where the experimental con- 
straints start to strongly limit the parameter space of mSUGRA. Given the upper limit a xp 
in mSUGRA, discovery of DM in the early stage of CDMS II experiment would enforce us 
to face to more flexible models . In this paper, we investigate the implications of recent 
measurement of muon anomalous magnetic moment for direct detection of Neutralino Dark 
Matter |32], |3(| |33|] in the various SUSY models. We impose the accelerator bounds of Higgs 
and sparticle masses. However, the prediction for b — > S7 decay rate is sensitive to details 
of the flavor structure of the soft breaking terms, unlike the quantities we consider here. We 
therefore do not attempt to analyze the constraint from b — > 37 decays quantitatively []. 

The muon anomalous magnetic moment will be further checked using four times more data 
collected already. On the other hand, it can be consistent with the standard model value if 
one takes into account other theoretical calculations of the hadronic vacuum polarization |35] . 

Hn mSUGRA, the constraint from b — > 37 decay rate is not important compared to that from Higgs mass 
for n > case Q . 
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In this regards, we consider two extreme cases for the range of Aa M ; 

27 x KT 10 < Aa M < 59 x 1(T 10 (la bound), (3) 

and 

< Aa M < 11 x 1(T 10 (2(7 below) (4) 

We study the minimal supergravity (mSUGRA) model as well as some other SUSY models 
where the assumption of strict scalar mass universality at the GUT scale is relaxed. Unlike 
mSUGRA, the relic density of neutralino could be very small for those models. Since the 
direct DM detection rate depends on the product of the a xp and local LSP density p x , one 
must consider what fraction of our local halo density p\ OC ai could be composed of neutralinos. 
Unless there is separation for different types of DM, the ratio of LSP DM to total DM should 
be the same locally in the Galaxy and globally in the whole Universe. In our paper we present 
the thermal relic density Q x and a xv separately. However, following assumption might be 
useful to estimate the signal ratio, 

P X = Plocal X (-^) (5) 
* 'DM 

where £Idm is the total contribution of DM to the total energy density of the Universe. When 
combined with the reduced Hubble constant h ~ 0.7 and VL DM ~ 0.3 we have VL DM h 2 ~ 0.15. 
For our present study, we assume that p x = pi oca i if ^\^ 2 — 0-15, and present the contours 
of future experimental reach using local LSP density rescaled as p x = pi oca i X (Q x h 2 /0.15) if 
Q x h 2 < 0.15. For convenience, we define the scaled cross section, 

oscaied = cr XP x (fi x /i 2 /0.15) for tt x h 2 < 0.15 

= a xp for tt x h 2 > 0.15. (6) 

We will find that the increase in the a xv may compensate the decrease in the Q x h 2 so that the 
DM detection rate more or less remains unchanged 0. For some model we also see that more 
than one order of magnitude enhancement of the counting ratio is possible when LSP is the 
subdominant component of the local halo. 



2 mSUGRA 

In the minimal supergravity model, it is usually assumed that all squared scalar masses receive 
a common soft SUSY breaking contribution m 2 Q = mfj = m 2 D = m 2 L = m 2 E = m 2 Hu = rn 2 Hd = 
m 2 at the GUT scale Mx — 2 • 10 16 GeV, while all gauginos receive a common mass M and all 
trilinear soft terms unify to A at the same scale. The renormalization group (RG) evolution 
of soft breaking squared Higgs masses then leads to consistent breaking of the electroweak 
symmetry, provided the higgsino mass parameter [i can be tuned independently [R7j. In this 
paper, we choose the weak scale input parameters mb(mb) = 4.2 GeV, m t (m t ) = 165 GeV, 
and tan (5. We minimize the tree level potential at renormalization scale Q = ^m^rn t , which 
essentially reproduces the correct value of p obtained by minimizing the full 1-loop effective 
potential [pq] . With these assumptions, the mSUGRA model allows four continuous free 



parameters (m, M, A and tan/?). We take p > because of BNL constraint. 
2 The model independent situation is recently discussed in |36|. 
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Figure 1: Allowed mass ranges of (a) chargino-sneutrino and (b) LSP-SU(2) singlet slepton in mSUGRA for 
tan/3 = 10. tan/3 = 50 case is shown in (c) and (d). We take A — and /i > 0, and scan m £ [0, 1] TeV and 
M £ [0, 1] TeV, subject to experimental constraints. The heavily marked points satisfy la bound of Aa M . 



For generic SUSY parameters, the chargino-sneutrino diagram provides a dominant contri- 
bution to a M . In Fig. 1 (a) and (c), we plot the allowed mass range of chargino and sneutrino 
for two different choices of tan/3 = 10 and 50 respectively. The allowed mass range of LSP 
and SU(2) singlet sleptons are shown in Fig. 1 (b) and (d) for tan/? = 10 and 50 respectively. 
Here, we take A = and ii > and allow m and M to vary in the intervals m < 1 TeV and 
M < 1 TeV. For the parameter range scanned in this plot, we find that the lightest Higgs 
boson h couples essentially like the single Higgs boson of the SM; we thus demand that its 
mass rrih > 111 GeV. This follows from recent LEP results [p9[] , allowing for a 2 GeV theo- 



retical uncertainty in the calculation of rrih |40|| . Further uncertainty on rrih comes from the 
error of top quark mass value (rrih is varied by ~ ±3 GeV when m t is varied by ±5 GeV 
around 175 GeV f4lf). We include loop corrections to the masses of neutral Higgs bosons 



from the third generation quarks and squarks, including leading two-loop corrections [40]. We 



further require that the chargino mass m~+ > 100 GeV |42|. We exclude regions where the 



LSP is charged, i.e., < m^o or < m^o. The parameter region is cut by Higgs mass 
constraint at left edge of the allowed region, and the lower edge of the allowed region is deter- 
mined by neutral LSP constraint. The heavily marked points satisfy the further requirement 
of 27 x 10~ 10 < Aa M < 59 x 10~ 10 . We use the expressions of ref.f28|j for the calculation of 

For tan/3 = 10, only small range of sparticle masses (m x ± ~ 200 GeV and ~ 200 GeV) 
are allowed by the lightest Higgs mass bound, the requirement of the neutral LSP and la 
bound of Aa M . For the small tan/? values, the lower bound of Aa M gives a strict upper limit 
on sparticle masses because Aa^ is proportional to tan/3 as mentioned already. On the other 
hand, the lightest Higgs mass bound rrih > HI GeV requires lower limit of sparticle masses 
because large radiative corrections from top and stop loops are required. These two opposite 
tendency allow only limited region of parameter space for the small tan/3 values. 
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Upper limit of sparticle masses increase as tan/3 increases. For tan/3 = 50, the maximal 
values of light chargino and sneutrino mass are ~ 430 GeV and ~ 820 GeV respectively within 
lcr bound of Aa M . For large tan/5 values, the region of parameter space where m x ±, irty is 
small is excluded by the upper bound of Aa M measurement because small sparticle masses 
gives too large Aa M values in this case. 

These upper (and lower) bounds of sparticle masses with the fact that \x > is preferred 
has a significant importance on the predictions of neutralino-proton cross section, o xv [|29], 
and neutralino dark matter relic density, Q x h 2 . In most situations the dominant contribution 
to the spin independent amplitude is the exchange of the two neutral CP-even Higgs bosons. 
The sizes of the XiXii^^H) couplings do depend quite significantly on the sign of /i, unless 
tan/? ^> 1 [I3|. In particular, for fi < strong cancellations occur [|30| both within different 
contributions to the same coupling, and between the h and H exchange contributions to a XP . 
Because positive /i is preferred by Aa M experiment, cancellations would not occur. Also the 
upper (lower) limit of sparticle masses has a significant effect on the minimum (maximum) 
values of neutralino-proton cross section and neutralino dark matter density. When the 1 a 




Figure 2: a xp vs. Q, x h 2 in mSUGRA (a) without and (b) with la bound of Aa M for tan/3 = 30 respectively. 
The corresponding plots are shown in (c) and (d) for tan/3 = 50 respectively. We take A = and \i > 0, and 
scan m € [0, 1] TeV and M £ [0, 1] TeV, subject to experimental constraints. The straight and dashed lines 
show a SC aied = 10 -8 and 10~ 9 pb, the approximate level of the claimed sensitivity by the CDMS (Soudan) and 
GENIUS proposals respectively. We also indicate the region of parameter space where Aa M is more than 2a 
lower than current central value ((b) and (d) lower region). 

bound of Aa^ is imposed, there is not too much parameter space remaining for tan /3 ~ 10 . In 
Fig. 2(a) ((c)) and (b)((d)), we therefore show the allowed region of a xp vs. Q x h 2 without and 
with Aa^ constraints for tan/3 = 30(50) respectively. Here, we take A = and \i > and allow 
m and M to vary in the intervals m < 1 TeV and M < 1 TeV and require rrih > 111 GeV, 
m x ± > 100 GeV and neutral LSP. The straight and dashed lines show the approximate level of 
the claimed sensitivity by the CDMS (Soudan) M and GENIUS |TD[ proposal respectively. The 
calculation of o xv is based on refs. [fMj, We use the value m s (p\ss\p) = 130 MeV for the 
strange quark's contribution to the nucleon mass; this matrix element is uncertain to about a 



5 



factor of 2, leading to a similar uncertainty in the prediction of a xp . Finally, the calculation 
of the scaled LSP relic density £l x h 2 uses results of refs.|43], |45|]; s— channel poles are treated 
as described in ref . [ 46] . The co-annihilation of x? with sleptons is not included. 



For tan/? = 30, a xp is larger than ~ 3 x 10~ 9 pb and can reach upto ~ 4 x 10~ 8 pb within 
la bound of Aa M (upper region in Fig. 2 (b)). In this case, Q x h 2 vary from ~ 0.2 to ~ 1.6. 
If we require a cosmological constraint Q x h 2 < 0.3 W7\, some portion of region are allowed 



and are within the approximate level of the claimed sensitivity by CDMS(Soundan) and the 
GENIUS proposal. Note that we do not include the fx? coannihilation here. The lower limit 
of Q x h 2 might be reduced substantially for the points of most left-down portion of parameter 
space in Fig. 2 (a) and (b) where parameter space is ended due to the condition that x? is n °t 
LSP. Note however that the prediction is very sensitive to the running of the third generation 



sparticle r between Planck scale and GUT scale [48, (49f|. For tan/? = 50, a xp is larger than 



~ 4 x 10~ 9 pb and can reach up to ~ 1 x 10~ 7 pb within lcr bound of Aa M . Q x h 2 is mostly 
less than 0.3 due to the effect of s-channel pseudoscalar Higgs pole f|3], [31], ^(J of neutralino 
pair annihilation and is thus cosmologically acceptable. The smaller value of itia is achived by 
the bottom Yukawa running of mn d enhanced by tan 2 (3. Heavy higgs exchange contribution 
in a xp is also enhanced as Hss coupling is enhanced and ran is reduced. Again all region is 
within the approximate level of the claimed sensitivity by the GENIUS proposal . 

On the other hand, when turns out to be more than 2a below the current central value, 
a XP is mostly below 1 x 10~ 9 pb and Q x h 2 is larger than 1(0.4) for tan/5 = 30(50). For 
tan/3 = 10, a XP < 0.5 x 10~ 8 pb and f2 x /i 2 > 0.4 is found in our numerical study. They are 
cosmologically disfavored unless fx? coannihilation is taken into account. 



3 More general models 

Now we relax our assumptions, allowing for non-universal soft scalar masses at the GUT scale, 
while keeping the unification of the gaugino masses. As specific models, we consider a model 



with non-universal Higgs mass and an SO (10) Grand Unified model |29 . 

In a model with non-universal Higgs, we assume that soft breaking Higgs mass terms are 
different from the other universal soft scalar mass terms at GUT scale, 

m ^ m Hu (= m Hd ) (7) 

For simplicity we keep m Hu = m Hd at the GUT scale []. 

The 5*0(10) theory incorporates a complete generation of MSSM matter superfields into 
the 16-dimensional spinor representation, \l/i6. In addition to these matter superfield, the 
minimal S'O(IO) model includes a 10-dimensional Higgs superfield $10 which contains the two 
Higgs superfields of the MSSM (as well as their SU(3) triplet, SU(2) singlet partners). When 
SO(10) breaks to the MSSM gauge group SU(3) C x SU(2) L x U{1) Y , additional D-term 
contributions (parameterized by M|> which can be either positive or negative) to the soft 
SUSY breaking masses arise pD[ : 

m Q = m l = m u = m w + M D 
m 2 D = m 2 L = m\ % - 3M 2 , m 2 Hu d = m 2 10 =F 2M 2 D , (8) 

where m\§ and mio are scalar soft breaking masses for fields in the 16 and 10 dimensional 
representations of S'O(IO), respectively. We consider this model too, but we assume mi 6 = 
mio = m for simplicity. 

3 A model with non-universal Higgs mass, eq. (0) can be considered as a special case of 50(10) model, eq. 
d) with M% = 0. 
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The modifications (0) and (]8|) of the mSUGRA boundary conditions change our predictions 
through the changed value of and pseudoscalar Higgs mass at the weak scale, which 
mostly affects o XP and Q x h 2 (and also through modifications of the slepton spectrum, which 
change Aa M and Q x h 2 , for a SO(10) model). The changes of value can be understood as 
follows. In the mSUGRA scenario, the contributions of m and M to the weak scale values of 
the soft breaking Higgs boson masses can be parameterized as 

m 2 Hd ~ m 2 + 0.5M 2 ; 
m 2 Hu ~ e H m 2 - 3.3M 2 , (9) 

where we have assumed sin j3 ~ 1 but ignored contributions from the bottom Yukawa coupling. 
The coefficient e# is small, because the GUT scale value of m 2 H is canceled by the scalar masses 
appearing in the RG running 0. The effect of non-universality on m 2 Hu can be parameterized 
by introducing roq + mjj + m 2 Hu = ?>m 2 and 5m 2 H = rn 2 Hu — m 2 : 

m 2 Hu =5m 2 H + e H m 2 s -3.3M 2 . (10) 

This equation follows because the radiative correction to ttVjj is proportional to m 2 q + mfj + 
m 2 Hu , hence the effect of RGE running is the same as in mSUGRA with the replacement 
m 2 — > m 2 . If sin (3 ~ 1, correct gauge symmetry breaking requires fi 2 ~ —m 2 H — M§/2, where 
all quantities are taken at the weak scale. In mSUGRA {i.e.,5m 2 H = 0), fi > M is predicted 
unless m is extremely large. However, if 5m 2 H is non-zero, it directly affects the value of at 
the weak scale. 

These modifications of the mSUGRA boundary conditions also have an impact on the 
Higgs masses. At tree level, the mass of the pseudoscalar Higgs boson is simply given by 

m A = m H u + m H d + 2 ^ 

~ m 2 Hd + » 2 -M 2 /2 (11) 

where all quantities are taken at the weak scale and in the second line, we used the relation 
fi 2 ~ ~ m2 H u ~ Mz/2 which holds if sin/? ~ 1. Therefore, m\ depends on rn 2 Hd and fi 2 at 
weak scale. When we increase the soft Higgs masses in same amounts, the increase of m 2 Hd 
is compensated by the decrease of fi 2 (through the increase of fn 2 Hu ). In the model with non- 
universal Higgs mass (with the assumption m 2 Hu = rn 2 Hd at GUT scale), m 2 A thus does not 
change too much along the change of m#/m. On the other hand, in the 5*0(10) model, the 
modifications of rn 2 Hu and m 2 Hd at GUT scale are in the opposite direction. In this model, 
m 2 A would be quite reduced by the decrease of m 2 Hd and fi when M|, has a large and negative 
value. 

The reduction of tra implies the reduction of charged Higgs mass m H ± and heaviest CP 
even Higgs mass m#. This has a significant implication on Q x h 2 , a xp and the branching ratio of 
b — > S7 decay. In mSUGRA model, LSP pair annihilation through XX ~ > A, H — > X channel is 
not important unless tan/3 3> 1. However, when rriA and m# decrease, the above channel can 
be quite enhanced especially in the resonance region. This in turn gives quite reduced Q x h 2 . 
On the other hand, the neutralino-proton cross section a xp increases as m# decreases. This is 
because the dominant contributions to a xv comes from the exchange of the two CP-even Higgs 
bosons. Also BR(b — > 57) is affected by m% value through charged Higgs-top quark loop. 

4 The effective eh is slightly negative for low SUSY breaking scale, but turns positive if this scale is large. 
Recall that the relevant scale for the analysis of gauge symmetry breaking increases with increasing sparticle 
masses. 
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The reduction of m% tends to increase BR(b — > 57) because the charged Higgs-top quark 
loop always interferes with SM contribution constructively. Therefore the preferred value of 
tan f3 is expected to be different from that of mSUGRA. Note that the ino-squark loop have 
a opposite sign to charged Higgs-top loop in mSUGRA model for fi > case. For very large 
tan/3, they tend to overshoot the charged Higgs contribution. 




Figure 3: a xp vs. Vt x h 2 in a model with non-unversal Higgs mass (a) without and (b) with Aa M constraint 
for tan/3 = 30 respectively. We take A = and fx > 0, and scan m E [0, 1] TeV , M S [0, 1] TeV and 
1 < m 2 Hu /m 2 < 20, subject to experimental constraints. The straight and dashed lines show the approximate 
level of the claimed sensitivity by the CDMS (Soudan) and GENIUS proposal respectively. 

Fig. 3 shows the correlation between a xp and Q x h 2 for tan/3 = 30, A = and fi > in 
a model with non-universal Higgs mass. Here, we scan the region m < 1 TeV, M < 1 TeV 
and 1 < m 2 Hu /m 2 < 20, under the same mass constraint in Fig. 2. In Fig. 3(a), the light 
(yellow) shaded region indicates m 2 Hu /m 2 = 1, i.e, mSUGRA case. The maximum value of a xp 
increases by about one order, compared to mSUGRA case and reaches up to ~ 7 x 10~ 7 pb 
where Q x h 2 value is about 0.02 and therefore LSP constitutes a subdominant part of local halo 
DM. The highest value of the scaled cross section, a sca ied is ~ 10~ 7 pb. In Fig. 3(b), we further 
require that the Aa M constraint is satisfied. In this figure, the dark (red) shaded regions are 
correspond to the la bound case and the medium (green) shaded regions to the more than 
2a below current central value case. For the la bound case, Q x h 2 values in large portion of 
allowed regions are less than 0.3 and most of such regions are within the approximate level 
of the claimed sensitivity by the CDMS (Soudan) proposal. Also, for the 2a below case, the 
large portion of allowed region are cosmologically acceptable and within the sensitivity level 
by the GENIUS, contrary to mSUGRA model where all regions are cosmologically disfavored 
and have very small a xp with 2a below case. 

Fig. 4 shows the correlation between a xp and Q x h 2 for tan/5 = 30, A = and fi > in 
5O(10) model. Here, we scan the region m < 1 TeV , M < 1 TeV and -50000GeV 2 < M^ < 
+50000GeV 2 , under the same mass constraint in Fig. 2 fj In Fig. 4(a), the light (yellow) 

5 We keep all soft scalar masses are positive at GUT scale, therefore — < M|, < ^- always. 
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Figure 4: a xv vs. il x h 2 in 50(10) model (a) without and (b) with Aa^ constraint for tan/3 = 30 respectively. 
We take A = and fi > 0, and scan m G [0, 1] TeV , M G [0, 1] TeV and -50000GeV 2 < < +50000GeV 2 , 
subject to experimental constraints. The corresponding plots are shown in (c) and (d) for fixed M5 = 
— 37500GeV 2 . We also show the scatter plot in (e) with Aa M constraint for tan/3 = 10. The straight and 
dashed lines show the approximate level of the claimed sensitivity by the CDMS (Soudan) and GENIUS proposal 
respectively. 
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shaded region indicates Mp = 0, i.e, mSUGRA case. The maximum value of a xp increases by 
more than two order of magnitude compared to mSUGRA case. It reaches up to ~ 2 x 10~ 5 
pb where Q x h 2 value is ~ 4 x 1CT 3 . The maximum value of cr S caied is ~ 2 x 1CT 6 pb f\ when 
VL x h 2 ~2x 10~ 2 . Therefore the LSP constitutes a subdominant part of the local halo DM in 
this case. In Fig. 4(b), we further require that Aa^ bound is satisfied. In this figure, the dark 
(red) shaded regions are correspond to the la bound case. Large portion of the parameter 
space have cosmologically allowed Q x h 2 < 0.3 and within the approximate level of the claimed 
sensitivity by the GENIUS proposal. 

For 50(10) model, the Higgs sector is changed drastically, as it can be read off directly 
from eq. (|3|). The heavy CP even and CP odd Higgs masses can be quite small if M|, has 
large negative values. In this CclSG , (J y^p £1X6 largely enhanced and the LSP sometimes hit the 
Higgs pole leading to greatly reduced Q x h 2 . There appear the region where the detection rate 
(which depends on the product o xv Q x ) are enhanced by more than one order of magnitude, 
compared to mSUGRA and non-universal Higgs mass model, within cosmologically acceptable 
region. This happens when <C 2m x , therefore pole effect is reduced while a xp is enhanced. 

Such highest detection rate comes from the points where LSP is a subdominant component 
of the halo Dark Matter, as can be seen in the figure. For clarity, we show the correlation 
between a xv and Q x h 2 in Fig. 4 (c) without and (d) with Aa M constraints for M|> = —37500 
GeV 2 . 

The medium (green) shaded regions correspond to the parameter space where a M is lower 
more than 2a compared to the current central value. Here the effect of the existence of the 
D term is small; a xp remains small and Q x h 2 is too high. This is due to our artificial cut 
off of the parameter space of M|,. To make /i small to have acceptable Q x h 2 , one needs M^, 
of roughly same order to that of gaugino masses. To compare we show the scatter plot for 
tan/? = 10 in Fig. 4 (e). The SUSY scale is substantially smaller compared to in Fig. 4 (b) 
(tan/3 = 30), therefore effect of M|, is visible. The la below case could have acceptable fl x h 2 
and a large counting ratio for any value of tan/5 when M|, is changed freely. 



4 Conclusions 

In this paper, we have investigated the implications of the recent measurement of muon MDM 
on the spin-independent neutralino-proton cross section a xv and LSP dark matter relic density 
fl x h 2 in several SUSY models with either universal or non-universal soft scalar masses at the 
GUT scale. We considered two extreme scenario for Aa M bound, i.e, 27 x 10 -10 < Aa^ < 
59 x 10~ 10 (la bound) and < Aa M < 11 x 10~ 10 (2a below). a xp and Aa^ can become 
large if tan/? is large, because both quantities are sensitive to chirality violation in the matter 
(s)fermion sector (which in the MSSM is enhanced for large tan/3) and if the SUSY mass scale 
is low. On the other hand, Q x h 2 become large if the SUSY mass scale is high. The la bound 
of Aa M prefer high tan/3 and low SUSY mass scale while for the 2a below case, low tan/3 and 
high SUSY mass scale are preferred. Therefore la bound of Aa M gives generally large a xp 
and small £l x h 2 while the opposite is true for 2a below. While these general statements are 
fairly model-independent as long as one sticks to the field content of the MSSM, quantitative 
predictions do depend significantly on details of the spectrum of superparticles, in particular 
on the implementation of radiative gauge symmetry breaking. 

We considered three different models: mSUGRA, a model with non-universal Higgs mass, 
and an 5*0(10) GUT model. In mSUGRA model, some portion of allowed parameter space 

6 This is the approximate level of claimed signal by DAMAp]. 



10 



is cosmologically acceptable. The a xp is within the sensitivity level by the GENIUS proposal, 
if we require la bound of Aa M . On the other hand, if Aa M is 2a below the current central 
value, the maximum a xp is mostly below 1CT 9 pb and Q x h 2 so large that it is cosmologically 
disfavored unless coannihilation of fxi is taken into account. 

In a model with non- universal Higgs mass (for 1 < m 2 Hu /m 2 and tan/3 = 30), the maximum 
value of a xp increases by about one order of magnitude while the minimum value of Q x h 2 is 
reduce by about one order of magnitude. With the la bound of Aa M , large portion of allowed 
parameter space gives cosmologically acceptable Q x h 2 and are within the sensitivity level of 
CDMS (Soudan) proposal. Even with the case where Aa M is 2a below from current central 
value, some portion of parameter space gives cosmologically acceptable Q x h 2 and are mostly 
within the sensitivity level of the GENIUS proposal. In this model, the maximal allowed value 
of dark matter detection rate (which depends on the product a XP Q x ) remains more or less 
same as that in mSUGRA model. 

Dramatic deviations from mSUGRA predictions are possible if one introduces 50(10) 
D— term contributions to scalar masses. In this model, within the la bound of Aa M , the 
maximal allowed value of a xp for given tan/3 can exceed the mSUGRA prediction by almost 
three order of magnitude, if the D— term contribution to the mass of the Higgs boson that 
couples to the top quark is positive and large. The corresponding Q x h 2 is also reduced, however 
the counting ratio maybe enhanced by factor of 10 compared to the previous models. In this 
model, the higgs sector is changed drastically. The heavy CP even and CP odd Higgs masses 
can be quite small if M|, has large negative values. The counting ratio is the highest when 
the LSP constitutes a subdominant potion of local halo DM. 

The predictions of a xp and Q x h 2 depend on the boundary conditions of model parameters 
at GUT scale. The FCNC constraints only require universality in squark sector in the same 
gauge multiplets. The universality of sfermions and higgs masses at the GUT scale as in 
mSUGRA requires /i 3> M, and it is not favored cosmologically. In mSUGRA model, only 
the parameter space with mass degeneracy of f and Xi or the 2m x ~ ttla is allowed. 

The too large relic density may imply that LSP might have more higgsino component than 
predicted in the mSUGRA model, or equivalently non-universal scalar masses at the GUT 
scale for higgs sectors. The improved measurement of Aa M , (non) discovery of dark matter 
at future search experiments, and improved observations of Q x h 2 will give us a hint on the 
LSP natures. Especially, if DM is discovered in the early stage of near future experiments, we 
should consider the variation of the soft mass terms. 

The implication to the collider physics is also important. Consistency to the cosmology 
force us to think about the parameter space where tan f3 ~ 50 or m x o ~ rrif in mSUGRA. For 
example, the degeneracy of stau and LSP mass indicates the domination of the decay mode 
into stau and very low acceptance of tau and anti-tau lepton pairs in x?, ~ > ^ T decays at LHC. 
The end point of m TT in the cascade decay is important for the determination of sparticle 
masses |||. However, it is equivalently possible that the nature's choice is a non-universal 
Higgs mass at GUT scale. In such case %° and xt niay have substantial wino component and 
therefore they will be frequently produced in q~L decay at LHC. The identification of such decay 
modes and the determination of fi parameter is demonstrated in |^TJ for the case where SUSY 
scale is relatively low and the cascade decay of ino to slepton has substantial branching ratios. 
When SUSY scale is high, the ino with substantial higgino and gaugino mixing dominantly 
decays into Higgs bosons. The study for such a case might be interesting to see the ability for 
LHC to measure the relation between fi and M to determine a xp and Q x h 2 . 
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